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Gravitating t'Hooft-Polyakov magnetic monopoles can be constructed when coupling the 
Georgi-Glashow model to gravitation. For a given value of the Higgs boson mass, these gravi- 
tating solitons exist up to a critical value of the ratio of the vector meson mass to the Planck 
mass. The critical solution is characterized by a degenerate horizon of the metric. As pointed 
out recently by Lue and Weinberg, two types of critical solutions can occur, depending on the 
value of the Higgs boson mass. Here we investigate this transition for dyons and show that 

OS ' the Lue and Weinberg phenomenon is favorized by the presence of the electric-charge degree of 

OS ' freedom. 

os : 

O ! I. INTRODUCTION 

CD . 

The Georgi-Glashow (GG) model consists of an SU(2) Yang-Mills field coupled to a (real) Higgs triplet. The 
C*~) ' self interation of the Higgs field by means of a Higgs potential leads to a spontaneous breakdown of the symmetry. 

One distinguished feature of this model is that it admits topological solitons : the celebrated t'Hooft-Polyakov 
r ^ j monopoles 

. By coupling the GG field theory to gravitation, one obtains the SU(2) Einstein-Georgi-Glashow model. Several 
£SJ ' years ago, it was shown that gravitating magnetic monopoles, as well as non-abelian black holes, exist in a certain 
CO ■ domain of the space of the physical parameters of the model . For instance for a fixed value of the Higgs boson 
mass, the gravitating monopoles exist up to a critical value a cl - of the parameter a (the ratio of the vector meson 
mass to the Planck mass). At the critical value a cr , a critical solution with a degenerate horizon is reached. In 
particular, for small values of the Higgs boson mass, the critical solution where a horizon first appears corresponds 
to an extremal Reissner-Nordstr0m (RN) solution outside the horizon while it is non-singular inside. 

Recently Lue and Weinberg reconsidered the equations of the self-gravitating magnetic monopoles and dis- 
covered an insofar not suspected phenomenon Indeed for large enough values of the Higgs boson mass, the 
critical solution is an extremal black hole with non-abelian hair and a mass less than the extremal RN value. An 
C"| , independent numerical analysis of the equations has confirmed the results of Lue and Weinberg and reinforced 
the agreement between the numerical solutions and the algebraic conditions these solutions have to obey. 
• i-H . During the last months, there was a growing interest devoted to monopole and dyons, considered as classical 
' solutions of the non-abelian Einstein-Born- Infeld-Higgs (EBIH) model ^ [u^. This is motivated by the fact that 
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the low energy effective action D-brane is related to a Born-Infeld model (see e.g. 11 1). Because Born-Infeld 
models naturally contain a theta-term proportional to (F^F^) 2 , their nonlinearities are better tested by the 
dyons solutions which excite both magnetic and electric fields. 

It is therefore natural to study if the gravitating dyons of the GG model are also sensitive to the Lue and 
Weinberg effect and to see if the degree of freedom confering the classical lump its electric charge favorizes or 
attenuates this effect. In this report we extend the analysis of Ref. [|| to the case of spherically symmetric 
gravitating dyons. We show that the transition also occurs for dyons and that the phenomenon is enhanced by 
the presence of the new field describing the electric charge degree of freedom. 



II. THE EQUATIONS 

We consider the SU(2) Einstein-Georgi-Glashow action In order to obtain static spherically symmetric 

globally regular solutions we employ Schwarzschild like coordinates for the metric 

ds 2 = -A 2 Ndt 2 + N^dr 2 + r 2 (d9 2 + sin 2 8d<p 2 ) , (1) 
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and introduce, as usual, the mass function m(r) by means of 

N(r)=l- 2 -^. (2) 

Then we use the standard Wu-Yang ansatz for the spatial components of the gauge field and the hedgehog ansatz 
for the Higgs fields (see e.g. j|Jq]) 

1 - Kir) - 1 - Kir) 

A r =0, A e = -e 4 , ^, A = eg ^sin6>, (3) 

9 9 

and 

$=e r H(r)v, (4) 

with the standard unit vectors e r , eg and while v denotes the Higgs field's vacuum expectation value. Since 
we want to describe dyon, we also employ the spherically symmetric ansatz for the time component of the gauge 
field 

Aq = vJ(r)e r . (5) 
It is also convenient to introduce the dimensionless coordinate x and mass function /i, 

x = gvr , /i = gvm , (6) 
as well as the dimensionless coupling constants a, f3 



a 



••-«** , 0-£, (7) 



where G is Newton's constant, Mfj is the Higgs boson mass and My/ is the gauge boson mass. 

With these ansatz and definitions, the classical equations of the model reduce to the following system of five 
coupled radial equations : jl^Jbi]] 



x 2j,2 j2 K 2 _ o 1 

/j' = a 



2 I , , NR >2 + _ Nx 2 H >2 



2 A 2 A 2 N 

+ {K l~2 1)2 + H2 K 2 + Y x2{H2 ~ 1)2 ) ' (8) 

and 

9 / 2J 2 K 2 2K' 2 „\ A 

for the functions parametrizing the geometry (the prime indicates the derivative with respect to x) ; for the matter 
functions we obtain the equations 

{ANK')i = AK(^-±+K--^\ , (10) 



x 2 J' 
~A~ 



2JK 2 
AN 



(11) 



and 
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{x z ANH')' = AH [2K 2 + ^-x 2 (H 2 -1)) . (12) 



The equations of motion depend on two physical parameters a, (3. This form of the parameters are used in |5J 
but several alternative notations were used in previous papers. In || they use a, b with 

a = 2a 2 , b = \{3 2 = X -f , (13) 

and in passing we redefined by /3 the parameter labelled (3 in |T^JT^| . 

In solving the equations, a, (3 are imposed by hand; the third physical parameter to be fixed by hand is the 
electric charge of the dyon, noted Q. It is encoded Jl^,[l5| into the asymptotic behaviour of the function J(x) 

J(x)\ x ^oo = Joo--+o{\) (14) 
X x z 

So that Q is imposed as a boundary condition of the function x 2 J 1 at x — oo. 

The regularity of the solution at the origin, the finiteness of the mass and the requirement that the metric (|l|) 
approaches the Minkowski metric for x — * oo lead to the following set of boundary conditions [Jill 

/i(0) = , K(0) = 1 , H(0) = , J(0) = (15) 



A(oo) = l , K(oc) = , ff(oo) = l , (x 2 J')U^oo = Q ■ (16) 

which fully specify the problem. 

In absence of gravity (i.e. a = 0, N = A = 1) the first two equations are trivial and the solutions are the dyons 
of Julia and Zee |fT2f . They were studied numerically in some details recently |fT5| . In particular, it was found 
that for 6 = dyons exist for an arbitrary value of Q while for b =/= the dyons exist only to a maximal value 
Q cr (see Fig. 8 of [|5|); for too high values of Q the value becomes larger than one and, as seen from Eq.©, 
K becomes oscillating. Within the range of (3 that we will explore in this paper (gravitating) dyons exist for 
Q < 0.7. 

The equations (|)-(|l|) were studied in details in Jl3|,[l4| for the case (3 — 0. In particular, it was shown that 
the non-Abelian gravitating dyon bifurcates at a — a cr into an extremal RN solution with 



N(x) = ^-*^+Q 2 ) 2 , A[x) = ! (17) 



K(x)=0 , H(x) = l , J(x) = Q(— =L_ - -) (18) 

ay/l + Q 2 x 



defined on x G [ct-y/l + Q 2 , oo]. In the next section we study the dyons solution for b > 0. 



III. GRAVITATING DYONS 



In order to make the paper self consistent, we first briefly recall how the magnetic monopole solutions approach 
critical solutions, when the vector boson mass, i.e. a, is varied, while the ratio of the Higgs boson mass to the 
vector boson mass, i.e. 6, is kept fixed. Recently, Lue and Weinberg Q realized that there are two regimes of b, 
each with its own type of critical solution. 

In the first regime b is small (b < 25), and the metric function N(r) of the monopole solutions possesses a single 
minimum. As the critical solution is approached, i.e. as a — > et cr , the minimum of the function N(r) decreases 
until it reaches zero at r — tq. The limiting solution corresponds to an extremal RN black hole solution with 
horizon radius = ro and unit magnetic charge for r > tq. Consequently, also the mass of the limiting solution 
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coincides with the mass of this extremal RN black hole. However, the limiting solution is not singular in the 
interior region, r < ro. We refer to this limiting approach as RN-type behaviour. 

In the second regime b is large, and the metric function N(r) of the monopolc solutions develops a second 
minimum as the critical solution is approached. This second minimum arises interior to the location of the first 
minimum, and decreases faster than the first minimum. Therefore, the critical solution is reached at r = r*, 
where the second minimum reaches zero. The critical solution thus possesses an extremal horizon at r* < ro, 
and corresponds to an extremal black hole with non-abelian hair and a mass less than the extremal RN value. 
Consequently, we refer to this second limiting approach as NA-type (non-abelian-type) behaviour. 

To summarize : the non-Abelian gravitating monopole exist on a portion of the (a, b) plane limited by a curve 
a cr (b). Somewhere on this curve (at (a cr = 3/2, b cr ~ 26.7)) there is a critical point separating the RN-type and 
NA-type of ending of the solution. 

By solving Eqs.(||) to (|lj) we were able to show that the same phenomenon occurs for dyon solutions. Before 
to present the details of the two regimes, we discuss the solutions in the case Q = 1/2, b = 50 which is generic of 
the NA-type (the corresponding critical a is a cr ss 0.83175). 

The profiles of the functions A, K, H at the approach of the critical value a cr was described at length in || 
and obey a similar pattern in the case of dyon. Therefore we present here the profiles of the functions N and J. 
We found it instructive to superpose on our figures the profiles of the functions 

i , B = — ~7=- (19) 
A AVN 



The ratio J/ A naturally appears when one eliminates the function A (by using (^)) from the other equations. The 
special combination B becomes a constant when the solution approaches an extremal RN black holes (p7|) , (|l8|) . 
Moreover it takes the value 

BRN = —7==S (20) 
ay/1 + Q 2 

this can be checked from the numerical solutions. 

So on Figs. 1-4, the functions N, J, J/ A, B are represented as functions of the variable y 

„=_-L_. (21) 
a^l + Q 2 

It is such that y a = 1 (with an obvious notation we also define y* as the value of the interior minimum of N) . 
Figure 1 is drawn for a = 0.75, i.e. much before a cr . The function N still present a unique minimum around 
y = 1. Figs. 2, 3 are obtained for a = 0.8315. Now N has developped a second minimum at y w 0.80. We 
observe that J (and then also J/ A and B) deviates only little from zero for y e [0, 0.8]. A little before y = 0.80 
the function B increases very quickly to attain the value (^pf ) around y = 1. We also notice that J/ A develops a 
ripple between the two minima on N. This general tendency of the functions is confirmed by Fig. 4 drawn very 
close to the critical a. The function B here increases very steeply (in fact quasi linearly) from zero at y* ~ 0.795 
to the value (|20| ) at y = 1. 

By solving the equations for several values of a,b and Q we obtain a strong numerical evidence that the two 
behaviours found in J(| persist in the case of dyons. This statement is illustrated by Figs. 5,6 where the ratio 
y* = r*/r characterizing the NA-type of behaviour is plotted respectively as a function Q (for a few fixed values 
of b) and as a function of of b (for Q = and Q — 1/2) [[l6|. The two figures indicates in particular that, when 
Q increases, the NA-type appears for lower values of b. 

Corresponding to Fig. 6, the critical value a cr is also plotted in Fig. 7. The remarkable fact indicated by this 
figure is that the line a cr (b) in the (a, b) plane is independent of Q (this was also observed in (m) for 6 = 0). Only 
the transition point (a cr , b cr ) is moving with Q on the line. 

The fact that the curve corresponding to Q = 1/2 overtakes the value a cr = 3/2 is not in contradiction with 
H (i.e. for Q — 0). Expanding the dyon equations around a double zero of the function N(r) (i.e. r — r* or 
r = ro) is the same way as in || would lead in principle to a Q-depending critical value of a limiting the RN 
and NA-types. Technically, it would involve the diagonalisation of a 3*3 matrix (rather than a 2*2) and leads to 
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a more involved algebra. However, owing that Q is related to the asymptotic decay of the function J it is very 
unlikely that an expansion of the solution around r* (or ro) could lead to an explicit expression for a cr (Q) and 
therefore we refrain from performing such an expansion. 

Clearly , Figs. 6, 7 report data which is obtained from a numerical analysis of the equation and probably the 
NA-type solution (with Q = or Q ^ 0) exists for still slightly lower values of b. These solutions cannot be 
constructed because of the limitation in accuracy of our numerical routine. In order to support this statement we 
report in Fig. 8 the value N m i n of the (RN-type) minimum of N when the internal minimum first appears, this 
as a function of b. As argued in [Q] we believe that the transition between RN and NA-types of behaviour occurs 
just when this curve crosses zero. An extrapolation of the curves gives b tr ~ 26.7 for Q = and b tr ~ 18.7 for 
Q = l/2. 



IV. CONCLUSION 



The Einstein-Georgi-Glashow model constitutes a good theoretical laboratory for testing the properties of 
gravitational solitons. The solutions are particularly rich and obey non trivial patterns of bifurcations of several 
types ]5|,[n3],[nj . Here we have analyzed numerically the equations for the spherically symmetric gravitating dyons 
in this model, focusing our attention on the intermediate region of the Higgs boson mass parameter. 

Our results gives strong evidence of two different scenarii of ending of the non-abelian dyons at some critical 
value a = a cr (b), confirming the recent result of |J. Three points desserve to be further stressed, (i) The domain 
of existence of gravitating dyons in the (a,b) plane seems to be independant of Q. (ii) The transition bewtween 
the RN and NA regime occurs at a Q-depending value (a C r,b cr ) of the parameters. For monopole (Q = 0) we 
have a cr = 3/2 , b cr « 26.7. Dyons have a cr > 3/2 and the precise value depends (rather indirectly) on the charge 
Q. (iii) At non-zero values of b the dyons solutions exist up to a maximal value Q C r(b), limiting the domain of 
investigation of the equations. 
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Figure captions 

Fig. 1. 

The profiles of the functions N, J, J/ A, B as functions of y for Q = 1/2, b = 50 and a = 0.75. 
Fig. 2. 

The profiles of the functions N, J, J/ A, B as functions of y for Q — 1/2, b — 50 and a = 0.8315. 
Fig. 3. 

Enlarged view of Fig. 2 in the region Y = 0.9 where the two minima of N occur. 
Fig. 4. 

The profiles of the functions N, J, J/ A, B as functions of y for Q = 1/2, b = 50 and a = 0.83175. 
Fig. 5 

The ratio = r»/r is presented as a function of Q in the NA-type regime for = 8 (b — 32) and for = 6.5 
(b = 21.125). 
Fig. 6 

The ratio y* = r*/ro is presented as a function of b in the NA-type regime for monopole solutions and for dyon 
solutions with charge Q = 1/2. 
Fig. 7 

The critical value a cr — 1 is presented as a function of b in the NA-type regime for monopole solutions and for 
dyon solutions with charge Q = 1/2. 
Fig. 8 

The value of the function N(r) at the first minimum where the second minimum first appears is presented as 
a function of the value of b, for monopole solutions and for dyon solutions with charge Q = 1/2. 
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Figure 4 
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